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The H ^ gg amplitude relevant for Higgs production via gluon fusion is computed in the 
four-dimensional helicity scheme (fdh) and in dimensional reduction (bred) at the two- 
loop level. The required renormalization is developed and described in detail, including 
the treatment of evanescent e-scalar contributions. In fdh and bred there are additional 
dimension-5 operators generating the Hgg vertices, where g can either be a gluon or an 
e-scalar. An appropriate operator basis is given and the operator mixing through renor¬ 
malization is described. The results of the present paper provide building blocks for further 
computations, and they allow to complete the study of the infrared divergence structure 
of two-loop amplitudes in FDH and bred. 
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1 Introduction 

Higgs production via gluon fusion is one of the most important LHC processes. Its com¬ 
putation at higher orders requires renormalization and factorization to cancel UV and IR 
divergences. The renormalization is less trivial than the one of standard QCD processes due 
to the required renormalization of non-renormalizable operators. The virtual corrections 
have been computed in conventional dimensional regularization (cdr) [1-5]; the required 
theory of operator renormalization in CDR has been developed in Ref. [6] , based on general 
work in Refs. [7, 8]. 

In the past years, several alternative regularization schemes have been developed. 
Purely four-dimensional schemes such as implicit regularization [9, 10] and fdr [11] have 
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been proposed and used to compute processes of practical interest such as H 77 [12, 13] 
and H ^ gg [14]. The present paper is devoted to regularization by dimensional reduction 
(dred) [15] and the related four-dimensional helicity (fdh) scheme [16]. Both schemes 
are actually the same regarding UV renormalization, but they differ in the treatment of 
external partons related to IR divergences.^ There has been significant progress in the un¬ 
derstanding of FDH and DRED: the equivalence to CDR [20, 21], mathematical consistency 
and the quantum action principle [22], infrared factorization [23, 24] have been established 
— these results solved several problems that had been reported earlier, related to violation 
of unitarity [25], Siegel’s inconsistency [26], and the factorization problem of [27, 28]. In 
addition, explicit multi-loop calculations have been carried out [29-33]. 

More recently, the multi-loop IR divergence structure of fdh and dred amplitudes has 
been studied in Ref. [34]. It has been shown that IR divergences in fdh and dred can be 
described by a generalization of the CDR formulas given in Refs. [35-38]. The description 
involves IR anomalous dimensions 7 ® for each parton type i. In Ref. [34] they have been 
computed for the cases of quarks and gluons by comparing the general IR factorization 
formulas with explicit results for the quark and gluon form factor. In FDH and DRED, 
however, the gluon can be decomposed into a T)-dimensional gluon g and {4 — D) additional 
degrees of freedom, so-called e-scalars g. In dred, e-scalars also appear as external states. 

The present paper is devoted to a detailed two-loop computation of the amplitude 
H ^ gg in FDH and DRED. In DRED, this involves the computations oi H ^ gg and 
H gg, since the external gluons can either be gauge helds or e-scalars. The fdh result 
is identical to the one for H ^ gg and has already been given in Ref. [34], but we will 
provide further details here. 

This detailed computation is of interest for two reasons: First, it provides the basis for 
obtaining the remaining IR anomalous dimension for e-scalars at the two-loop level. Second, 
it provides an example of the required renormalization in fdh and dred, including operator 
renormalization and operator mixing. The difficulty of renormalization in FDH and DRED, 
particularly in connection with H —> gg, has been pointed out e. g. in Refs. [33, 39]. 

The outline of the paper is as follows: Section 2 gives a brief description of the regular¬ 
ization schemes and of the relevant Lagrangian and operators. It ends with a detailed list 
of the required ingredients of the calculation. Apart from the actual two-loop computation 
and ordinary parameter and field renormalization that are described in Sections 3 and 4, 
respectively, the main difficulty lies in the renormalization and mixing of the operators 
generating H —)• gg. This is discussed in general in Section 5, and specific two-loop results 
are presented in Section 6 . Section 7 then provides the final results for the on-shell ampli¬ 
tudes for H ^ gg and H —>■ gg. The appendix contains details on our projection operators 
and gives Feynman rules for the different operator insertions. 

®Parts of the literature, e. g. Refs. [17-19] used the term DR/dimensional reduction for what is called 
FDH here. 
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CDR 

HV 

FDH 

DRED 

internal gluon 

gbU 

gbU 

gbU 

gb- 

external gluon 

gb^ 

-gbU 

-gbl' 

gbU 


Table 1. Treatment of internal and external gluons in the four different regularization schemes, 
i.e. prescription which metric tensor has to be used in propagator numerators and polarization 
sums. 

2 Regularization schemes and H ^ gg 

It is useful to distinguish the following regularization schemes [24]: conventional dimen¬ 
sional regularization (cdr), the’t Hooft-Veltman (hv) scheme, the four-dimensional helic- 
ity (fdh) scheme, and dimensional reduction (dred). In all these schemes, momenta are 
treated in D = 4 —2e dimensions (the associated space is denoted by QDS with metric ten¬ 
sor In order to dehne the schemes, one also needs an additional quasi-4-dimensional 

space {Q4:S, metric g^'^) and the original 4-dimensional space (dS, metric g^'^). The treat¬ 
ment of gluons in the four schemes is given in Tab. 1. In the table, “internal” gluons are 
defined as either virtual gluons that are part of a one-particle irreducible loop diagram or, 
for real correction diagrams, gluons in the initial or final state that are collinear or soft. 
“External gluons” are defined as all other gluons. 

Mathematical consistency and D-dimensional gauge invariance require that QAS D 
QDS D 45" and forbid to identify g^'^ and g^'^. Details can be found in Refs. [22, 24, 34], 
The most important relations for the present paper are 

+ = = = ( 2 . 1 ) 

where a complementary 2e-dimensional metric tensor has been introduced. With these 
metric tensors we can decompose a quasi-4-dimensional gluon field as 

At^ = + gf^^A^ = > + (2.2) 

into a D-dimensional gauge field A^ and an associated e-scalar field A^ with multiplic¬ 
ity W = 2e. ^ Correspondingly, there are two types of particles in the regularized theory: 
D-dimensional gluons g and e-scalars g. The unregularized external gluons g of fdh are a 
part of g. 

The regularized Lagrangian of massless QCD then reads 

^QCD, regularized = “ ^{9^ A^,af + i + d^^CaD^Ca + (2.3a) 

A = —{b^A-ubA.,). - qAm - 4 A^AfiAAs.d- 

(2,3b) 

^In many applications of fdh the dimensionality of QAS is left as a variable Ds, which is eventually set 
to Ds = 4. The multiplicity of e-scalars is then W = Ds — D. 
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Here, and = d^+iggA^ denote the non-abelian field strength tensor and the covari¬ 
ant derivative in D dimensions; "0 and c are the quark and ghost fields. In Eq. (2.3b) the 
coupling of e-scalars to (anti-)quarks is given by the evanescent Yukawa-like coupling g^. 
This could in principle be set equal to the strong coupling gg. But, since both couplings 
renormalize differently this would only hold at tree-level and for one particular renormal¬ 
ization scale [20]; the same is true for the quartic e-scalar coupling g^^. In Eq. (2.3b) 
we introduce an abbreviation that includes the appearing Lorentz and color structure; 
{die)abed sleifa.befcded'^'*9^^ -|-perm.), where “perm.” denotes the 5 permutations aris¬ 
ing from symmetrization in the multi-indices {a, a) ... (0,7). In the following we use all 

„2 

couplings in the form cti = 47 with i = s, e, 4e. 

The process H ^ gg is generated by an effective Lagrangian which arises from in¬ 
tegrating out the top quark in the Standard Model. In CDR it contains only the term 
— ^XHFa'^F^y^a- In fdh and dred one again has to distinguish several gauge invariant 
structures containing either D-dimensional gluons or e-scalars. The effective Lagrangian 
can be written as 


Teff = XHOi + X^HOi + Y, (2.4) 

i 

with 

Oi = (2.5a) 

Oi = -l(D^i-),(D^i,),. (2.5b) 

denote operators involving products of four e-scalars. Such operators are not impor¬ 
tant in the present paper and will not be given explicitly. Like for ag^oie and the 
couplings A and can be set equal at tree-level, but they renormalize differently and have 
different (3 functions. 

Our final goal is the calculation of the two-loop form factors for gluons and e-scalars. 
This requires the on-shell calculation of the 3-point function T hai^A'' {<1^ —p, —r). All mo¬ 
menta are defined as incoming, so q = p + r. The 3-point function can be separated 
into and corresponding to the amplitudes for H ^ gg and H -A gg, 

respectively. In dred, both on-shell amplitudes are needed according to Tab. 1. In fdh, 
only H ^ gg is needed, which however is identical io H ^ gg and will not be discussed 
seperately. 

The on-shell calculation requires the knowledge of the two-loop renormalization con¬ 
stants and These in turn can be obtained from an off-shell calculation of 

rffAi^A‘'- Projectors extracting the required renormalization constants from the off-shell 
Green functions and precisely defining the gluon and e-scalar form factors are given in 
appendix A.l. 

We have now all ingredients to discuss the classes of Eeynman diagrams that contribute 
to Fj/ai^A'' ill FDH and dred: 
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Figure 1. Sample two-loop diagrams for the process H ^ gg and H ^ gg in dred. e-scalars are 
denoted by dashed lines. The appearing coupling combinations from left to right are Xa^, X^a^, 

1. Genuine two-loop diagrams Some remarks concerning the calculation are 

presented in Sec. 3. 

2. Counterterm diagrams arising from one- and two-loop renor¬ 

malization of the helds, the gauge parameter and of the couplings a^, Og, and 04 ^. 
The required renormalization constants are presented in Sec. 4. 

3. Counterterm diagrams ^-rising from one-loop renormalization of the effective 

Lagrangian (2.4) at the one-loop level, which includes the renormalization of A and 
A^. This is a major complication and will be presented in Sec. 5. 

2LCX b 

4 . Overall two-loop counterterm diagrams ^haiiX’' a-rising from the two-loop renormal¬ 
ization of the effective Lagrangian (2.4), equivalently from the renormalization con¬ 
stants and 5Z'^. These renormalization constants are generally defined by the 
requirement that the appropriate off-shell Green functions are UV finite after renor¬ 
malization. For the case of an elegant alternative determination is possible [6], 
but that method fails for 5Z\^. The results for SZ"^ and iZ"^ are presented in Sec. 6. 

3 Genuine two-loop diagrams 

As mentioned above the Green function can be separated into and ^ 

corresponding to H ^ gg and H -A- gg. Examples for genuine two-loop diagrams with 
either external gluons or e-scalars are shown in Fig. 1. 

All loop calculations have been performed using the following setup: the genera¬ 
tion of diagrams and analytical expressions is done with the Mathematica package Feyn- 
Arts [40]; to cope with the extended Lorentz structure in QiS we use a modified version 
of TRACER [41]; all planar on-shell integrals are reduced and evaluated with an inple- 
mentation of an in-house algorithm that is based on integration-by-parts methods and the 
Laporta-algorithm [42]; all non-planar and off-shell integrals are reduced and evaluated 
with the packages EIRE [43] and FIESTA [44]. 

4 Parameter and field renormalization in FDH and DRED 

We now consider the counterterm contributions and They are given by 

diagrams exemplified in Fig. 2, where the counterterm insertions are generated by the usual 
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H 


H 



Figure 2. Sample one-loop counterterm diagrams originating from the renormalization of the 
couplings as,ae, a 4 e, and of the gauge parameter respectively. 


multiplicative QCD renormalization of the couplings and fields present in Eq. (2.3b). In 
the following we present the values of the required (3 functions and anomalous dimensions, 
which govern the renormalization constants. 


4.1 (3 functions 

The renormalization of the couplings , Oe > and 04 ^ is done by replacing the bare couplings 
with the renormalized ones. As renormalization scheme we choose a modified version of 
the MS scheme: like in Ref. [34] we treat the multiplicity of the e-scalars as an initially 
arbitrary quantity and subtract divergences of the form ■ As a consequence, the 

corresponding (3 functions depend on N^: (^) = a 4 e 5 A^e), with 

i = s,e, 4e. They are given in Refs. [33, 34] and read: 




—CaNf - 2CfNf 


(4.1a) 


Os \ 2 / Uf 


47r / V 47r 


CfNfN, 


+ C>(q;^), 




Ca(2 - N,) + Cf(-4 + N,) - Nf 


+ 0{a^). (4.1b) 


The renormalization of the quartic coupling {<^4e)abcd more complicated since the tree- 
level color structure, fate fade, is not preserved under renormalization [20]. In the case of 
an SU(3) gauge group one therefore has to introduce three quartic couplings, with 
i = 1,2,3, each of them related to one specific color structure in a basis of color space. 
Examples for such a basis are given e. g. in Refs. [29, 30]. 

In the present case oi H ^ gg the renormalization constant for only appears in 
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diagrams like the third of Fig. 2. Hence, only the following contracted /3 function is needed: 


Jabcd 


^ab Safj 



+ 



2 


CaNf{ 4: -2N,) + CrNpi-S 


4iV,) 


+ (1 - - 2iV.)| + 0{a^). 


(4.2) 


This result is obtained from a direct off-shell calculation. It agrees with a general result 
from [45]. 


4.2 Anomalous dimensions 


For the off-shell calculation of Thai^A'' also renormalization of the fields and of the gauge 
parameter ^ is needed. The renormalization of ^ is fixed by the requirement that the gauge 
fixing term does not renormalize: ^ —>■ The anomalous dimensions 7 ^ = /r^^^lnZj of 

gluon and e-scalar fields are obtained from a direct off-shell calculation of the respective 
two-loop self energies. Their values up to two-loop level read: 




47r/ 


dvr/ 


'6 2 6 / 3 






- -^CaNf - 2CfNf 


(4.3a) 


7d 


Ois 

47r 


)(7a(3-0 


Qe 

dvr 


dvr/ 


dvr/ 


C 


61 


- Nf 


11 


-a'"'’"- 


CaNf ( —1 + + CfNf ^2 -|- — 


CXs 

dvr 


CTe 

dvr 


— 5 Cf^f 




+ 0{a^). 


(4.3b) 


Setting Nf: and Oe to zero in Eq. (4.3a) yields the well-known gluon anomalous dimension 
in CDR, see e.g. [46]. The value of 7 ^ agrees with the general result for the anomalous 
dimension of a scalar field [45], confirming the point of view that e-scalars behave like 
ordinary scalar helds with multiplicity 
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5 Operator renormalization and mixing in FDH and DRED 

The second type of counterterm contributions, denoted by originates 

from the necessary renormalization of the effective Lagrangian (2.4), equivalently of the 
operators Oi and Oi. One major difficulty is that multiplicative renormalization of the 
parameters A and is not sufficient since the operators mix with further operators. We 
will show that the full operator mixing involving gauge non-invariant operators has to 
be taken into account. The renormalization constants cannot be predicted from known 
QCD renormalization constants but need to be determined from an off-shell calculation. 
The general theory of operator mixing in gauge theories and the classification of gauge 
invariant and gauge non-invariant operators has been developed long ago [7, 8 , 47]. 

In the following we briefly describe operator mixing in the much simpler case of CDR 
and then explain the cases of fdh and dred, which involve further operators. 


5.1 Operators in CDR 

In CDR, a useful basis of scalar dimension-4 operators, which is closed under renormaliza¬ 


tion, is given in Ref. [6]: 

Oi = (5.1a) 

02=0, (5.1b) 

O3 = ^i’0i’, (5.1c) 

O 4 = - {d^Ca){d^Ca), (5.Id) 

O 5 = {D^^dfjC)aCa. (5.1e) 


Operator Oi is gauge invariant and related to coupling renormalization; O 2 = rmp'tp in 
Ref. [ 6 ] and corresponds to the fermion mass renormalization; we set m = 0. All other 
operators are constrained by BRS invariance and Slavnov-Taylor identities [7, 8 ]; operators 
O 4 and O 5 are not gauge invariant. The basis is chosen such that O 3 , O 4 and O 5 are related 
to field renormalization of '0, and c, respectively. In particular, the first two terms of 
O 4 are generated by applying the functional derivative 


A:ix) 


5 

<5i^(x) 


(5.2) 


on the gauge invariant part of the QCD action; the remaining term is then required by 
BRS invariance and the non-renormalization of the gauge fixing term.^ 

®See Refs. [ 8 , 47] for more details; the full operator O4 can be obtained from evaluating + 

W(d‘'ca)A,^^a, where W is the linearized Slavnov-Taylor operator and is the source of the BRS trans¬ 
formation of AA^ in the functional integral. Since W is nilpotent, this definition shows that O 4 is compatible 
with BRS invariance and the Slavnov-Taylor identity and can appear in the operator mixing. 


9 




The operators renormalize as 


Oi y ^ijOj^\,Qxei (^•^) 

where Oj^bare arises from Oj by replacing all parameters and fields by the respective bare 
quantities. Following an elegant proof in Ref. [ 6 ] the nontrivial CDR renormalization matrix 
Zij can be written in the form 

Zij = 6 ij + Di InZj. (5.4) 

Here, JD* are derivatives with respect to parameters and Zj are combinations of ordinary 
QCD renormalization constants. As a result, in particular the renormalization of Zn is 
given by 

Q 

Zii = l + as-—InZa^, (5.5) 

UOtg 

with the multiplicative renormalization constant of , Za^. In this way the renormalization 
of the parameter A in the CDR version of Teff is related to the renormalization of as- 

5.2 Operators in FDH and DRED 

In FDH and DRED, the basis of operators needs to contain additional terms involving e- 
scalars. We use a basis constructed analogously to Eqs. (5.1) from gauge invariant operators 
and operators corresponding to field renormalization. Then there are two kinds of changes: 
there are modifications of the operators O 3 and O 4 , and there are additional basis elements. 
The new basis operators correspond to the e-scalar kinetic term, Oi, to the new parameters 
cte and 040 O 3 and 04 e,j, and to the field renormalization of A^, O 4 . The notation is chosen 


such that in all cases Oj and Oj have a similar structure: 

Oi = (5.6a) 

O 2 = 0, (5.6b) 

03 = ^ip I/) ip- (5.6c) 

O 4 = Al{b^^F^y)a + gsfabc{d^K)-^li,bAy^c “ “ {9^Co) (^^Ca) , (5.6d) 

O 5 = {b^d^Ca)Ca, (5.6e) 

61 = (5.6f) 

O 3 = gei^M, (5.6g) 

04 = i^(.D^.D^i.)a, (5.6h) 

64 ,,^ = 0{Ab. (5.6i) 
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Since we consider massless QCD there is no e-scalar mass term. Like in Eq. (2.4), operators 
involving four e-scalars are not needed explicitly. 

This set of operators differs in a crucial way from the CDR case. The difference be¬ 
tween operators Oi and O 4 is related to the total derivative OA^A^^. Hence, the basis for 
space-time integrated operators (zero-momentum insertions) does not coincide with the 
one for non-integrated operators (non-vanishing momentum insertions). As discussed by 
Spiridonov in Ref. [6], in such a case his method cannot be used. Therefore, in fdh and 
DRED it is not possible to derive complete results for the operator mixing analogous to 
Eqs. (5.4) and (5.5). 

This implies two difficulties: Eirst, the two-loop renormalization of Oi and the corre¬ 
sponding parameter cannot be obtained from a priori known two-loop QCD renormal¬ 
ization constants but need to be determined from an explicit two-loop off-shell calculation. 
Second, the off-shell Green functions get contributions from unphysical, gauge non-invariant 
operators, so the full operator mixing needs to be taken into account. 

We have carried out the explicit one-loop calculations to obtain all required one-loop 
results for and Z^j. The results are 



(5.7a) 

(5.7b) 

(5.7c) 

(5.7d) 

(5.7e) 

(5.7f) 

(5.7g) 

(5.7h) 

(5.71) 
(5.7j) 
(5.7k) 

(5.71) 


Renormalization constants involving operators O 3 or are not needed for the calcula¬ 
tions in the present paper. The renormalization constants (5.7a)-(5.7d) agree with those 
given in Ref. [34]. The only gauge-dependent quantity is Z^g. This is due to the fact that 
operator O 4 is related to the field renormalization of the e-scalars. In all other renormal- 
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Figure 3. Sample one-loop counterterm diagrams originating from operators O3, O 4 , O 4 and O5. 


ization constants related to field renormalization the gauge-dependent parts incidentally 
cancel out. 

With these results the bare effective Lagrangian can be written as 

^ E bare + A, bare) , (5-8) 

3 

where the sum runs over all operators in Eqs. (5.6). Sometimes it is useful to write this 
using multiplicative renormalization constants for A and as 

= ^AAi/Ol,bare + A.i/Oi^^are + ' ' ' > (5-9) 

suppressing operators not present at tree level, such that XZx = XZu + X^Zi^ and similar 
for Za,. 

The one-loop counterterm effective Lagrangian involving the renormalization constants 
of Eqs. (5.7) is then given by 

= hY,[x dZlj^Oj + A, SZ-l^O^ ) . (5.10) 

3 

We have now all ingredients for the one-loop counterterm diagrams relevant for 

the computation of 77 —?■ gg, where the gluons are either H-dimensional gauge fields or 
e-scalars. These counterterm contributions arise from one-loop counterterm diagrams with 
one insertion of Sample diagrams are given in Eig. 3. They show insertions of oper¬ 

ators O 3 , O 4 , O 4 and O 5 . Feynman rules for operator insertions are given in appendix A. 2. 

The calculation shows that all these operators generate non-vanishing contributions 
to However, in the extraction of the form factors and two-loop renormalization 

constants to be discussed in the next section there are cancellations, and O 4 is the only 
new operator which contributes. 
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6 Two-loop renormalization constants of A and A, 


Putting together the results from the previous three sections it is possible to calculate 
the two-loop renormalization constants and appearing in Eq. (5.9). They can 

be obtained from a complete off-shell two-loop calculation and the requirement that the 
corresponding Green-functions are UV finite after renormalization: 


-p2L I -plLCT.a , p2LCT,a plLCT,b p2LCT,b 

J- ^ HAt^A’^ + -*■ HAi^A'' + HAi^A'' + ^ HAt^A'' 


off-shell 


= 0 . 

UV div. 


( 6 . 1 ) 


All ingredients except the last term are computed in the previous sections, and Eq. (6.1) 
is then used to extract and The result for ^Z"^ is: 


5Z 


2 L^/^\ 


dvr/ 




m _ n AT , N‘. 

Q Q “T 


36 _|_ 3 ^ 


-|- CaNf 


-# + 'iN, 


10 


+ CfNf- + 


( 6 . 2 ) 




Since the off-shell calculations have been done numerically with the help of FIESTA [44] 
the analytical expressions have been obtained by rounding to a least common denominator. 
The numerical uncertainty is less than ^ for the terms of the order 0(e~^) and | for the 
terms of the order 0 (e~^). 

Result (6.2) is not new; it agrees with Ref. [34], where it has been obtained using 
Spiridonov’s method. The recalculation serves as a test of the setup and the results given 
in the previous sections. At the same time a comparison with Ref. [34] confirms that 
Eq. (6.2) is actually exactly correct, in spite of numerical uncertainties. 

In the same way, we obtain the renormalization constant 5Z^: 


5Z' 


2 L _ / “s y 




f+ fiV. , + + 

o T 


A 


N, 


-\- Ca^f 


- 2 + 


a.s \ f a. 




q 3 

_ £ 2 +’^A, 

2 ' 

e 


CfNf 


3 i-3^ 


- 2 + 




6 , 

o “r 


\caNf 


_3 I 3jy 

O' A-‘-^e 




-|- CfNf 


3N, 3-In, 

2e2 ^ e 




iV2 


_ A A 




-l-JV. , 15 

e2 ^ 8e 


1 5 _ 2 — 

1 6 ^Ae 


(6.3) 
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Compared to Eq. (6.2) this result is more complicated and includes all combinations of 
the three couplings as^a^ and a/^^. This result is new; as described in Sec. 5 it cannot 
be obtained using Spiridonov’s method. The numerical uncertainty is less than ^ for all 
terms. A forthcoming comparison with a prediction of the infrared structure ol H ^ gg 
will conhrm that expression (6.3) is exactly correct [48]. 


7 UV renormalized form factors of gluons and e-scalars 

Now that all renormalization constants are known it is possible to calculate the two-loop 
form factors of gluons and e-scalars in the fdh and dred scheme. We present the results 
in two ways: First, we give results with independent couplings needed to determine the IR 
anomalous dimensions of gluons and e-scalars; second, we give simplified results, where all 
couplings are set equal. These can be viewed as the final results for the UV renormalized 
but IR regularized form factors. We give them including higher orders in the e-expansion. 


7.1 Results for independent couplings 

The UV renormalized but IR divergent form factor for H ^ gg ui DRED is given at the 
one-loop and two-loop level by 






11 _L V 
3 6 


e 




2 + H^(3) + 3^w) 




2 Nf 


+ 0 {e^), 


(7.1) 


F?^{as,ae,Xe/^,N^) 




-T + 


6 




175 
18 6 


1 + 2 ^ + 


N't 

36 


_ llvr2 _ Mtf3) + iV ( 49 , ZR _ 29 
27 36'' 3 ^ 27 ^ 72 3 




49 TTf _ 

27 72 3 A 


^ 6 A 


-|- CaNf 


7 -13 2.AT M Zl!. 2 Ae /Y 

' _|_ 3 ^ _|_ 27 ^ 18 4^ 3 A 


3e3 


+ CfNfI + ^ 


_ / a 

Vdvr 




(7.2) 


As mentioned in the beginning the g form factor in DRED is identical to the gluon form 
factor in FDH, and Eq. (7.2) agrees with the result given in Ref. [34]. 

Since there are no external e-scalars in diagrams related to the gluon form factor 
internal e-scalars have to be part of a closed e-scalar loop or have to couple to a closed 
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fermion loop. Hence, the effective coupling always appears together with at least one 
power of in Eqs. (7.1) and (7.2). 

The e-scalar form factor for H ^ gg in DRED is given by 


(Q^s, CTe, Cr4e, A/Ag, A^g) 



vr 


--T -2 + - + 2 - + e( - 4 + Y^ +yC(3)+4- 


TT 
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(7.4) 


Compared to Eqs. (7.1) and (7.2) the result with external e-scalars is more complicated and 
includes all combinations of the couplings as,ae and a 4 e. In this result, like in all previous 
results, the evanescent coupling ae appears always together with at least one power of Np 
and the quartic coupling a 4 e is always accompanied by a factor (1 — A^)- 

7.2 Results for equal couplings 

During the renormalization process the couplings as, ae, a/^e and A, A^ have to be distin¬ 
guished. After renormalization they can be set equal, giving a simpler form of the hnal 
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result.^ The results for = 2e at the one(two)-loop level up to order 0{e^) (0(e^)) then 
read: 
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"^If the results of Sec. 7.1 were not desired for independent couplings, the genuine two-loop diagrams 
could have been computed in a simpler way, with all couplings set equal from the beginning — this is what 
is done in many applications of fdh and dred in the literature. 
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8 Conclusions 


We have computed the H ^ gg amplitudes at the two-loop level in the fdh and dred 
scheme and presented the MS renormalized on-shell results up to the order e^. In dred, 
this involves two different amplitudes for H ^ gg and H ^ gg with external gluons/e- 
scalars. The computation is motivated because it contains key elements which constitute 
important building blocks for further computations, and because it is essential for the 
complete understanding of the infrared divergence structure of edh and dred amplitudes. 

The renormalization procedure has been described in detail. It is less trivial than in 
many QCD calculations in CDR, since not only the strong coupling needs to be renormalized 
but also evanescent couplings of the e-scalar. The computation provides a further example 
of the well-known fact that regardless of whether fdh or dred is used, the evanescent 
couplings have to be renormalized independently. 

Further, the renormalization of the effective dimension-5 operators involves mixing 
with new, e-scalar dependent operators. A suitable basis of operators has been provided. 
One unavoidable fact is that the extended operator space contains operators which are total 
derivatives. As a result the required operator mixing renormalization constants cannot be 
obtained in the same elegant way of Ref. [6] as in CDR. Instead, they had to be obtained 
from explicit one- and two-loop off-shell calculations. 

The results for the UV renormalized but infrared divergent form factors can also be 
used to complete the study of the general infrared divergence structure of two-loop ampli¬ 
tudes in FDH and dred, begun in Ref. [33, 34]. From general principles it is known that all 
infrared divergences can be expressed in terms of cusp and parton anomalous dimensions. 
The results of the present paper allow to extract the final missing two-loop anomalous 
dimension for external e-scalars. This extraction, together with further checks and results, 
will be presented in a forthcoming paper [48], where the infrared structure will also be 
investigated by a SCET approach. 
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A Appendix 

A.l Projectors and form factors of gluons and e-scalars 

According to its Lorentz structure the on-shell Green-function can be represented 

as 

=a{p-r) + bp'^rf^ + cp^^r'' + dp>^p’^ + (A.l) 

where the coefficients o... e are momentum-dependent quantities, and coefficient a is the 
gluon form factor. Due to QCD Ward-identities the relation a = —b holds, see e. g. Ref. [1]. 
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Accordingly, the on-shell Green-function with external e-scalars can be represented 

as 


-'on-shell 


= fip-r) 


(A.2) 


HAi^A^ 

where we refer to / as e-scalar form factor. All coefficients of the covariant decomposition 
can be extracted with appropriate projection operators that are given below. 

In the off-shell case the UV divergence structure of represented in a 

more specific way as 
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{p ■ r) + B p^'r^ -|- Cp^r’’' + Dp^p^ -|- E r^r^', 
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where the coefficients A... E are now momentum-independent. Since these divergences can 
be absorbed by counterterms corresponding to operators Oi and O 4 the relation A = —B 
again holds, see e. g. Feynman rules (A. 6 ) and (A. 12 ). Due to this there are two possibilities 
of extracting coefficient A, which corresponds to the desired renormalization constant 
The first one is to extract the coefficient of {p ■ r) g^'^ and neglect terms oc the second 

is to extract coefficient —B. We checked explicitly that the relations a = —b and A = —B 
hold throughout the paper. 

Again, the covariant decomposition with external e-scalars is much simpler and reads: 
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(A.4) 


The desired coefficient for the computation of is E. Accordingly, we extract the 

coefficient of (p • r) g^'^ and neglect terms oc p^, r^. 

The corresponding projection operators are: 

- {p'^rf^ p^r''){p ■ r) (A.5a) 
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(A.5c) 


A.2 Feynman rnles 

In the following we give Feynman rules according to operators Oi, Oi, O 4 and O 4 that are 
needed for the renormalization in the fdh and dred scheme. Feynman rules including four 
e-scalars are not relevant in this paper and are not given explicitly. 
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Feynman rules according to the Lagrangian term \HOi: 
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• Feynman rules according to the Lagrangian term Xf^HOi 
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Feynman rules according to the Lagrangian term H0^\ 
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• Feynman rules according to the Lagrangian term H0^\ 
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